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Abstract 

We study the ability to control d-dimensional quantum systems (qudits) encoded in the hyper- 
fine spin of alkali-metal atoms through the application of radio- and microwave-frequency magnetic 
fields in the presence of inhomogeneities in amplitude and detuning. Such a capability is essential 
to the design of robust pulses that mitigate the effects of experimental uncertainty and also for 
application to tomographic addressing of particular members of an extended ensemble. We study 
the problem of preparing an arbitrary state in the Hilbert space from an initial fiducial state. We 
prove that inhomogeneous control of qudit ensembles is possible based on a semi-analytic protocol 
that synthesizes the target through a sequence of alternating rf and microwave-driven SU(2) ro- 
tations in overlapping irreducible subspaces. Several examples of robust control are studied, and 
the semi-analytic protocol is compared to a brute force, full numerical search. For small inhomo- 
geneities, < 1%, both approaches achieve average fidelities greater than 0.99, but the brute force 
approach performs superiorly, reaching high fidelities in shorter times and capable of handling in- 
homogeneities well beyond experimental uncertainty. The full numerical search is also applied to 
tomographic addressing whereby two different nonclassical states of the spin are produced in two 
halves of the ensemble. 
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I. INTRODUCTION 



Motivated by applications ranging from controlling chemical dynamics to quantum infor- 
mation processing (QIP), control over quantum systems has become an increasingly impor- 
tant tool [THU]. In the context of QIP, the majority of studies of such processes are based on 
a collection of two-level (d = 2) subsystems or qubits, a natural extension of binary classical 
logic. In practice, the physical objects that encode quantum information never have solely 
two levels and control is necessary to isolate a particular qubit of interest. Alternatively, 
higher dimensional d > 2 systems, or qudits, can be employed for base-d quantum logic. 
While from a computer science perspective encoding in qudits does not lead to a change 
in computation complexity, from a physical perspective, information processing with qudits 
provides different trade-offs. For example, a log 2 d reduction in the number of subsystems 
must be compared with the complexity in implementing the universal SU(d) gates and 
quantum error correcting codes [7j. Moreover, qudits are interesting in their own right ex- 
hibiting properties such as nonlocality without entanglement [8] and providing a platform 
for explorations of quantum chaos [9]. 

A natural qudit is the hyperfme manifold of magnetic sublevels associated with the 
ground-electronic state of atoms, providing a Hilbert space of dimension d = (2J+ 1)(2/+ 1) 
where J is the electron angular momentum and I is the nuclear spin. Motivated by ongoing 
experiments in the Jessen group [TU], we work with 133 Cs, an alkali-metal atom with one 
valence electron J = S — 1/2 and nuclear spin / = 7/2, yielding hyperfme coupled spins 
of magnitude F = I ± S = 3,4, and a total Hilbert space of dimension d = 16. Exten- 
sions to other elements, such as the rare-earths, open opportunities to even larger Hilbert 
spaces (e.g., and d = 128 in holmium pT|). Early proposals for qudit control in hyperfme 
manifolds involved the use of a series of two-photon Raman transitions between magnetic 
sublevels [12] . More recently, time-dependent magnetic fields and a static tensor light shift 
have been used for arbitrary state preparation in the lower hyperfme manifold [10J. Alter- 
native schemes that employ radio-frequency and microwave magnetic fields have also been 
studied for applications in qudit quantum control P [131 E]- 

Our goal in this article is to expand the control tool box for hyperfme qudits, with 
attention to alkali-metal atoms, the standard elements used in laser cooling experiments. 
We particularly develop the methods of ensemble control, previously studied for spin- 1/2 

3 



nuclei in the context of NMR |15j . whereby subsystems are subjected to an ensemble of 
different Hamiltonians. Such methods are important in a number of scenarios. In the context 
of experimental uncertainty in the control parameters, the different dynamics generated 
by different members of the ensemble correspond to errors. Making the desired dynamics 
insensitive to this uncertainty is known as "robust control," and will be essential for the 
levels of precision required in implementing QIP protocols. In addition, an inhomogeneous 
ensemble may be intentionally imposed. In such a scenario different members of an ensemble 
may be tomographically addressed in space or time. 

The idea of robust control goes back to the development of composite pulse techniques for 
spin- 1/2 nuclei in NMR [161 HZ! where special sequences such as "CORPSE" and "SCRO- 
FULOUS" were designed to perform particular SU(2) rotations in a manner that is robust 
to errors in detuning and/or Rabi frequency. In more recent studies, Khaneja and Glaser 
showed that one can achieve control of spin-1/2 particles as a nearly arbitrary function of 
detuning or Rabi frequency [151 [181 EE]. Researchers are also applying such NMR-inspired 
tools to the control of cold atoms. For example, the groups of Blatt and Chuang used 
robust pulses to improve the performance of quantum algorithms with trapped ions [20J, 
and Jessen's group explored robust control neutral atom qubits trapped in optical lattices 
[2*T] . Tomographic addressing of ultracold atoms is becoming increasingly important. Using 
real magnetic fields, the Bloch group addressed the "shells" of a Mott insulator of |22j and 
the Meschede group addressed sites of a one dimension lattice lattice beyond the diffraction 
limit [23]. Even higher resolution addressing is possible using the spin-dependent light shifts 
of off-resonant laser fields |24j . Using the light shift of a focused laser beam to create a large 
spatial gradient, Bloch's group individually addressed atoms in a 2D lattice with a spacing 
closer than the diffraction limit of the addressing beam [25J. Such tomographic addressing 
can benefit from more sophisticated control analysis. 

In this article we study robust control and tomographic addressing of hyperfine qudits 
within the unified framework of ensemble control. In Sec. Ill] we describe the physical system 



and the available control Hamiltonians, and use these in Sec. Ill A as the basis for a semi- 
analytic protocol for synthesizing arbitrary states through a series of SU(2) rotations on 
overlapping subspaces. Based on the known results in inhomogeneous control of SU(2) 
rotations, we extend our qudit state preparation routine to the case of ensembles. With the 
semi-analytics in hand, in Sec. |III B we present an alternative approach to state synthesis 
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FIG. 1. Radio- frequency (rf) magnetic fields (in red) and microwave magnetic fields (in blue) 
control the ground-electronic hyperfine manifold of magnetic sublevels in 133 Cs. SU(2) rotations 
in the irreducible subspaces, F = 3, 4, are driven by rf-magnetic fields. Because of the sign difference 
of the (/-factors for the upper and lower manifold, the two rotations are in opposite directions. The 
microwave field is shown to be resonant only with the transition between the states \F = 4, mp = 4) 
and \F = 3,niF = 3) due to the application of a bias magnetic field that breaks the degeneracy 
between the magnetic sublevels. Microwave-driven SU(2) rotations on this qubit, together with the 
rf-control, can generate an arbitrary unitary transformation on the full 16-dimensional manifold. 



based on numerical optimization. We compare the two approaches by exploring robust state 
preparation in the presence of inhomogeneities. Beyond robust control, we study how we can 
employ the tools of ensemble control to spatially address different regions of a cold atomic 
cloud and perform local state preparation to create highly nonclassical, nonequilibrium states 
of the gas. In Sec. |V]we summarize and present the outlook for future research. 
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II. CONTROL HAMILTONIAN 



We will focus on control of the spin state of an alkali-metal atom in its electronic ground 
state through magnetic interactions. The governing Hamiltonian of the system is 

H = AI-S + (g s fi B S- gi fx N I)-B(t), (1) 

where the first term represents the hyperfine coupling between the atom's nuclear spin, I, 
and the valence electron's spin, S, and the second represents the interaction of the spins 
with the controlling magnetic fields. We consider three contributions to the fields, 

B(t) = B e s + B rf {t)+B ltw (t), (2) 

whose effect on the energy levels is depicted in Fig. ([I]). The first contribution is a static 
bias magnetic field which splits the energies of the magnetic sublevels in the linear Zeeman 
regime, HbBq <C A, while the next two terms are control fields oscillating at radio- (rf) and 
microwave-frequencies that drive transitions between those levels. We work in a regime in 
which the hyperfine coupling is significantly stronger than the interaction due to the applied 
magnetic fields. Then it is convenient to split the state space into a direct sum of spaces 
with total angular momentum f^> — I ± |. In the linear Zeeman regime, we can use the 
Lande projection theorem to write its contribution as 

Hb = Yl ^bBq-F^. (3) 

/=+,- 

Defining the Zeeman frequency f2 = —Q-^bBq (here and throughout % — 1), the total 
static Hamiltonian becomes 

H = ^p(P« - PH) + n (g r F^ - Fj")), (4) 

where AEhf is the hyperfine splitting, and are the projectors on the hyperfine manifolds 
f^>. To take into account the small difference in magnitude and opposite signs of the g- 
factors in the lower and upper manifold arising from the nuclear magneton, we have defined 
9r — \g+/9-\- The rf-field resonantly couples magnetic sublevels within a subspace with a 
given total angular momentum. Defining the Larmor frequencies Q x (y) = —2g-^BB x ( y ), we 
can again use the Lande projection theorem to write the rf-Hamiltonian as 

H rf = 2tt x cos(u r f t - <j) x )(g r Fi +) - F^) + 2Vl y cos(oo rf t - - F^'). (5) 
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Transforming into a rotating frame according to U r f = e l0J rft( F z f z )^ anc j ma ki n g the 
rotating wave approximation, yields 

H r f + H =g r (Q x cos <\> x - Q y sin <j) y ) Fj + g r (Q x sin <f) x + Q y cos <j) y ) 

— (Q x cos0 x + Q y sin0y) F^ + (Q x sin(j) x — Q y cos<f) y ) F^ 

A W 

+ A^pM - Fj~)) + (1 - *) Wr/ F« + ^(PW - PW), (6) 

where A = Q — w r /. 

The microwave power, frequency, and polarization are chosen so as to couple only two 
magnetic sublevels between the hyperfine manifolds, leaving the other states unaffected due 
to off-resonance effects. These states define a pseudo-spin- 1/2 |Jf) = \ f <y+ \m^ = , ||) = 
\f(~',Tn^ = f^)} with the usual Pauli operators, and the microwave-driven Hamiltonian 
takes the standard spin-resonance form, 

Hyw = COSfu^t - 4>nw)^x- (7) 

We now transform into a total rotating frame given by Utotai = £/ r /e~"**( p(+)_p( > - )/ ' 2 with 
a = u^ w — {m^+m^)uj r f and make the rotating wave approximation. Choosing the microwave 
frequency to be resonant with the pseudo-spin transition, = (g+rn^ — g^m^iiBBo + 
Ehf, the total Hamiltonian in the rotating frame, including the static rf and microwave 
contributions is 

H = g r [Q x cos 4> x - Vt y sin 4> y ] F^ +) + g r [tt x sin <p x + Q y cos <p y ] F y (+) 
— (Q x cos<fi x + Q y s'm(f) y ) Pj - ' + (Q x s'm(f) x — Q y cos</> y ) F^> 
+ (cos <p^ w cr x - sin^o-j,) - ]-(g r - l)m t u rf F^ +) (P M - P (-) ). (8) 

The control parameters available in this system are the amplitudes and phases of the two rf- 
and microwave-fields, {fl x ,(j) x ,fl y ,(j) y , 0,^,(1)^}. With these controls, we have shown that 
one can synthesize any unitary transformation on the hyperfine manifold [6l [13] . 

If we choose the phases of the x and y rf-coils to be such that <fr x = <p y — tt/2 = <p r f and 
choose the powers in the two coils to be equal, Q x = Q y = Q = Q r f, then the rf polarization 
is positive-helicity circular and the rf-field is resonant only with the lower hyperfine manifold, 
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leaving the upper manifold fixed. The Hamiltonian then reduces to 

H = H rf + H uw (9a) 
H rf = 2tt rf (- cos (p rf F^ + sin <j> rf F y H) , (9b) 
H uw = (cos (j)^ a x - sin <p^ w a y ) . (9c) 

Through our choice of amplitudes and phases of the applied fields, we have restricted 
the dynamics to an 8D subspace spanned by the 7D basis in the F = 3 manifold plus 
a single auxiliary state in the upper manifold, \F = 4, rrtp = 4) = The independent 
rf or microwave Hamiltonians, H r f and H uw , are each generators of an SU(2) rotation in 
an irreducible subspace. As we will show in the next section, with the given controls, a 
unitary transformation on the 8D Hilbert space can be constructed from a sequence of 
SU (2) rotations. 



III. CONTROL PROTOCOLS 

A. Semi- Analytic State Synthesis 

We study the problem of synthesizing an arbitrary state within the 8D Hilbert space 
defined above through a series of SU(2) rotations generated by the control Hamiltonian, 
Eq. Q. As our fiducial initial state, we begin with all the population in As a proof 
of principle, we employ the technique originally developed by Law and Eberly [26] in which 
we solve the inverse problem - begin with an arbitrary target state \ipT) in the space and 
then map it to the fidulcial state, If each of the controls can be reversed, we can then 
determine the mapping |^) — > \i[>t)- We will see that unlike the original Eberly and Law 
protocol that involves a qubit, this method will only perform approximate state mapping, 
with an error that decreases exponentially with the length of the pulse. 

To find a map that transfers the target state I^t) - ► It); we solve a sequence of max- 
imization problems. First we find an rf-pulse that maximizes the population in the state 
\F = 3,rriF = 3) = |!) and then we find a microwave-driven rotation that maximizes the 
population of \F = 4,rriF = 4) = The latter transformation can be found analytically 
because every unitary transformation in this 2D subspace is an SU(2) rotation. A state 

|X) = ct It) + cj, II) (10) 
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is defined by a Bloch vector, r, for the pseudo-spin. If n is a vector that bisects the z-axis 
of the Bloch sphere and r, then a 7r-rotation around n will drive all the population in \x) to 
|t). For the rf-driven rotation, the situation is slightly more complex. The 577(2) matrices 
represent only a subgroup of the general SU(2F + 1) unitary transformations on the space. 
Define an SU(2) spin coherent state for F — 3, parametrized by a direction on the sphere, 
\0,4>) = K^{6,(f)) ||). We we thus seek the F = 3 irrep of the SU(2) rotation TZ^(9, <p) 
that satisfies 



This is equivalent to the maximum value of the Husimi distribution with respect to spin- 
coherent states, Q(9,<fi) = \(9, 4>\4>t)\ 2 i which we can find easily. 

How optimal is this procedure in the achievable fidelity and the required number of steps? 
Since the microwave-driven rotation can completely transfer all population from the \\) state 
to |t), the only question is how much population in the lower manifold we can transfer to \\) 
using rf-driven rotations. Because the Husimi distribution Q(6, <fi) is everywhere positive, 
we can find a lower bound by looking at the case where the Husimi distribution is flat, i.e., 
the maximally mixed state. For a spin-F the Husmi function in such a case would have a 
uniform height of 1/{2F +1). Therefore the amount of population remaining in the lower 
manifold, and thus our error, is bounded from below by (2F/(2F + l)) n where n is the 
number of iterations. This implies that the population decreases exponentially with n and 
we can rapidly map \ipr) — > It)- By reversing the sequence of inverse-rotations, we achieve 
the desired state preparation, |t) — > \^t)- 

Because this state synthesis protocol is based on a series of SU(2) rotations, we can 
draw on previously known results in SU(2) control to extend our construction to the case 
of inhomogeneous parameters. In particular, Kobzar et al. [TH] used a powerful numer- 
ical routine to find pulse sequences that synthesized arbitrary unitary transformations in 
two- level systems as a function of inhomogeneous controls, and Li and Khaneja explored 
the theoretical underpinnings |18j. The key result is that for two-level systems driven by 
canonical Hamiltonians of the form 



max | (||ft (3) (fl,0) W) 



(11) 



H[e, A, Q(t), (j){t)] = (1 + e)Q(t) (cos</>(t) <r x + sin <f)(t) a y ) + Act 



(12) 
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one can synthesize arbitrary 577(2) transformation, 




) 



(13) 



V Jo 



as arbitrary functions of e > —1 and A. Here Q(t) and 4>(t) is are the amplitude and phase 
of the drive waveform that act to control the system. Whereas we allow for inhomogeneity 
in the amplitude, the phase is assumed to be known and controlled precisely. This result 
generalizes trivially for any irreducible representation of SU(2) on a spin of magnitude F. 

This construction implies that in principle one can design robust waveforms so that 
U (e, A) is essentially flat, and that one can design tomographic waveforms that address 
different members of the ensemble for well-chosen values of e and/or A. In practice, the 
controllability of the system will be limited by practical constraints such as bandwidth, slew 
rates, and time of interaction. The duration of the waveform is of particular importance 
given the ultimate constraint of decoherence. Note that the waveform that generates a 
desired U(e,A), is not unique. Whereas the proof given by Li and Khaneja provides an 
analytic algorithm to synthesize U(e,A), in practice, numerical optimization schemes lead 
to faster pulse sequences while maintaining very high-fidelity. 

We study here an extension of the results of ensemble control of SU(2) rotations to our 
protocol for qudit control. The evolution of each member of the ensemble is described by 
the extension of Hamiltonian, Eq. rt9J), 



Our goal is to synthesize different target states as a function of e MTO , e r f, and A. For each 
SU(2) rotation we can draw on the results of Li and Khaneja. We must ensure that the 
effect of inhomogeneity still allows efficient coherent transfer of probability amplitude from 
F — 3 to and that controllability is respected in the sequence of rf- and microwave-driven 
rotations. 

Because our protocol interleaves pulses acting on different subspaces, we must carefully 
examine the controllablity with respect to different inhomogeneities. Consider first the case 
of ensemble control in the presence of microwave amplitude inhomogeneity alone. The first 




(14) 
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rf-pulse must maximize the population in \\), averaged over e^ w , so Eq. (11) becomes 



max^KM^Me^))! 2 . (15) 
u,<p — 

€ \X'W 

Because the rf pulses cannot distinguish between the different members of this particular 
ensemble of states, this is equivalent to maximizing the population in \\) given an initial 
state, Peff, where 

Peff = N^2\i/>r(e ftw ))(il) T (e im )\, (16) 

with N as a normalization constant. In the worst case p e ff is a completely mixed state. In 
this situation the rf pulse accomplishes nothing, and 1/7 of the population remains in |4,). 
The first microwave 7r-pulse will transfers all of the population in |i) to The second rf 
pulse acts on another incoherent mixture over e^ w , but not a completely mixed state because 
we have emptied all of the population in In the worst case, 1/7 of the remaining 6/7 
of the population will be transferred by the rf-pulse to \\). The second microwave pulse in 



the sequence will in general act on an ensemble of spinors of the form in Eq. (10), \x( e ^w)}, 
because there now exits an ensemble of complex amplitudes in the superposition between 

and |t). However, because we can synthesize an arbitrary SU(2) matrix as a function 
of e^, we can find a transformation that rotates \x{ e ^w)) to |^) for the relevant range of 
values of e^ w . By repeating this procedure, we continually increase the population in \\) at 
an exponential rate, even in the presence of microwave power inhomogeneity. 

The case of rf inhomogeneities alone, \ipT{^rf)), is less favorable. Based on the results 
of Li and Khaneja [18], we can find an SU(2) rotation which maximizes the population in 
\l) for all e r f, and the first microwave 7r-pulse will transfer this population to |t). Whereas 
the application of the second rf pulse will again maximize the remaining population in 
F = 3 to the state \\) for all e r /, the second microwave pulse will generally not work 
as needed. Because the microwave field cannot distinguish between different spinors in 
Eq. (10) as a function of e r f, |x(e r /)), the state seen by the microwaves is mixed, p e ff = 



NJ2 e f \x( e rf))(x( e rf)\- in the worst case, this could be either a state whose Bloch vector 
points along the z-axis of the Bloch sphere or is the completely mixed state. In either case, 
the microwaves cannot, on average, increase the amount of population in the \f) state. Thus, 
without further care, this procedure cannot synthesis arbitrary states for different e r /. For 
the particular case that the target state is independent of e r f, we can overcome this difficulty 
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since we can synthesize a target state at the end of each rf pulse that satisfies this property. 
By employing such robust control pulses, we ensure a state \x) of the pseudospin that can 
be rotated by microwaves to |t) at each stage. But, more general cases are problematic. 

Finally, we consider inhomogeneities in the detuning as might arise from an external 
magnetic field (either noise or intentionally applied). In this case we will need to consider 
the role of reversibility in the state synthesis routine with more care. According to the Law 
and Eberly protocol, we design a pulse sequence that maps the target \iJjt) to the fiducial 
state according to ||) = Un---U 2 Ui \^t)- Then the desired state preparation follows as 
\4>t) = Ull/l-.-Ulf |t)- The inverse unitary is generated by the Hamiltonian — H. We can 
achieve this by simply adding 7r to the phase control waveform <fi(t) and also inverting the 
detuning, A — > —A. Thus, if we want to synthesize a target |^r(A)) we must use a sequence 
of rotations that maps \iPt(— A)) — > |t) to search for the control waveforms. 

In the presence of detuning inhomogeneity, there is an additional source of errors due to 
phases accumulated in the rotating frame. For instance, when the rf pulse is driving the lower 
manifold for time t r f, the population in the |t) state will acquire a phase = g r Am^t r f. 
We must compensate for this phase in the design of subsequent microwave pulses. This 
is possible since we can synthesize arbitrary microwave pulses as a function of detuning. 
Similarly, we can compensate for A-dependent phases accumulated in the F = 3 manifold 
while the microwave pulses are applied through the design of subsequent rf-pulses. 

Similar arguments may be made when multiple inhomogeneities are present. Based on 
these arguments we find that arbitrary target states of the form \i^t(^hw, A)) may be synthe- 
sized. We emphasize that although the particular protocol presented here cannot synthesize 
states that are different for different e r f, this does not imply that such pulse sequences do 
not exist. The existence of a semi-analytic protocol capable of synthesizing states which 
vary with e r / is a matter of ongoing research. 

B. Fully- numerical state synthesis 

Although the results of the previous section indicate that it is possible to perform qudit 
ensemble control with a sequence of well-designed SU (2) rotations, the duration of the pulse 
sequence can be unnecessarily long. We therefore also consider fully numerical optimization 
to search for the desired controls. The semi-analytic protocol lends confidence that such 
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a search routine will yield high fidelity results. For fully numerical searches we allow the 
microwaves and rf fields to be applied simultaneously and thus the evolution is no longer 
a simple series of SU(2) rotations on the given subspaces. We have found that we can 
attain full control of the system using only piecewise constant variations in the phase. This 
has the advantage that we may perform the search via unconstrained optimizations that 
often converge more rapidly than constrained routines. Searches then return a vector of the 
sequence of control parameters over N steps, 

= (0r/,l, (j>fj,w,l, 4>rf,2, 0/m>,2, 4>rf,N, 4>^w,n)- (17) 

To perform the numerical search for control pulses, we choose as our objective function the 
fidelity 

J-($)= l(^T(e^,A)|[/(e r/ ,e^,A,$)|t)| 2 , (18) 

where [^(e^, A)) is the target state for a given value of the ensemble. Unconstrained 
optimization is performed using the Matlab routine fminunc. 



IV. RESULTS 



A. Robust control: semi-analytical vs. fully numerical search 

We present results from several examples of robust control so that the target state is 



independent of e^, e r f } and A. The Hamiltonian is described by Eq. (14) with Q, /2ti = 
100 kHz. The maximum rf-Larmor frequency is VL x ^/2tx = 1.5 kHz, and the maximum 
microwave-Rabi frequency is Q fiw /27r = 3.5 kHz. For these parameters, we respect the 
linear Zeeman regime and can safely ignore off-resonant effects of the microwave radiation, 
in line with our model of dynamics on an 8D Hilbert space - the span of the F = 3 subspace 
plus |t). Given these time scales, we take 125 /is steps for piecewise constant evolution of 
the control waveforms. 

We compare the performance of the semi-analytical approach of Sec. Ill with the fully 
numerical search. In all cases, we performed optimizations for 20 states chosen according to 
the Haar of measure on SU(8) and then averaged the fidelity over the results for different 
states. In the semi-analytical approach, we alternate SU(2) rf- and microwave-generated 
rotations. Each rotation is made robust by decomposing them into 3 steps of 125 /is duration 
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whose amplitude and phases are found by numerical optimization. As the protocol does not 
have a specified end time, the number of rotations is not fixed. For each state we repeat the 
optimization 10 times and pick the optimization that attains a chosen threshold fidelity in 
the shortest time. In contrast, for the fully numerical approach, the microwave and rf fields 
are present simultaneously and we fix the total duration of the waveform, which sets the 
number of 125 /is steps. We choose the duration to be long enough to achieve high fidelities; 
the duration is chosen for each case depending on the amount of inhomogeneity. We set the 
amplitude of the rf and microwave fields to their maximum values, and only optimize their 
phase. The algorithm is then iterated until a target fidelity of 0.99 is reached. 

A comparison of the results for errors of up to 1% in e r /, e^ w and A are shown in Fig. [2j 
In both the semi-analytical and fully numerical protocols, the pulse sequences are found 
by optimizing the fidelity on a coarse grid in parameter space defined by e r j = 0, ±0.01, 
e fj, w = 0, ±0.01 and A = 0, ±0.01f2 r /, and then averaged. To ensure that the pulse sequences 
perform as desired, we then calculate the fidelity on a finer grid with 15 evenly spaced points 
between ±1% in those parameters. For the semi-analytical case, the total time to reach a 
fidelity greater than 0.99 is 3.94 ms, averaged over the 20 states. As can be seen from the 
figure, a fidelity of over 0.99 is maintained over the range of parameters, and the fidelity 
averaged over that range is 0.997. For the fully numerical case, as seen in Fig. [2j a fidelity of 
over 0.99 is maintained throughout the parameter range, and the fidelity averaged over that 
range is 0.994. While this sample shows slightly lower average fidelity, the duration of the 
fully numerical waveform is only 1 ms long, nearly a factor of 4 speed up compared to the 
sequence of SU(2) rotations. Thus, we see that for inhomogeneities of 1%, both approaches 
achieve high fidelity, though the sequence found via fully numerical optimization can do so 
in significantly less time. 

While the semi-analytical method yield waveforms that require substantially more time 
to achieve high fidelity when compared with the fully-numerical method, it serves as an 
important proof-of-principle. For the remainder of the paper we focus on the fully-numeric 
method. Figure ([3ja) shows the results of fully numerical optimization with errors of 5%. 
In this case, the optimization was performed on a grid in parameter space defined by e r / = 
0,±0.05, e^ w = 0,±0.05 and A = 0,±0.05fi r /. The duration of the entire pulse sequence 
is 5 ms. After the optimization, we calculate the fidelity on a finer grid of 15 points evenly 
spaced between ±5% for each of those errors, as shown in the figures. A fidelity of over 
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(a) 



(b) 

1 




FIG. 2. Average fidelity vs. inhomogeneity for (a) semi-analytic state preparation and (b) fully 
numerical state preparation. The x-axis represents the fractional variation in the each parameter 
for the three cases - power inhomogeneity in rf and micorwave, e r fn W , and detuning as a fraction 
of 0,nw We show the performance as a function of one inhomogeneity and averaged over the the 
other two. Each curve is an average of 20 random states in the 8 dimensional Hilbert space. In the 
semi-analytic case, each waveform consists of a series of alternating rf and microwave SU{2) pulses 
for a total of 3.98 ms. The optimized pulses are found through a numerical search of the time 
varying amplitude and phase of the relevant fields. In the fully numerical case, rf and microwave 
fields are applied simultaneously for a total of 1 ms. The microwave and rf powers are set to 
their maximum values, and a we perform a numerical optimization of the phases. For both the 
semi-analytic and fully numerical case, the waveform consisted of a series of piecewise constant 125 
Us pulses. 



0.988 is maintained throughout this range of parameters and we achieve an average fidelity 
of 0.993. Optimization by the fully numerical method is sufficiently powerful to compensate 
errors well beyond experimental uncertainty. As an example, in Fig. (J3jb) we present the 
results of optimizing errors of A that are 10% of Q r f. We perform the optimization on a grid 
of 10 evenly space points between ±0.1Q r f. The timing is the same as the previous example. 
Fidelities of over 0.98 are maintained over the entire range of the detuning, while the fidelity 
averaged over the full range of errors is 0.992. We summarize our results in Table [1} 
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FIG. 3. (a) Fidelity of fully numerical state preparation vs 5% errors in e r j, e^ w , and A/Q r f. 
The x-axis represents the fractional variation in the each parameter for the three cases - power 
inhomogeneity in rf and microwave, tvf^wi and detuning as a fraction of ^l^ w . We show the 
performance as a function of one inhomogeneity and averaged over the the other two. Each curve 
is an average of 20 random states in the 8 dimensional Hilbert space, (b) Fidelity of fully numerical 
state preparation vs 10% errors in A/£l r f. 



Inhomogeneity Range 


Optimization Method 


Total Time 


Average Fidelity 


e rf ,e^ w ,A/n rf £ [-.01, .01] 


Semi- Analytic 


3.94 ms 


0.997 


e rf ,e^ w ,A/n rf G [-.01, .01] 


Fully Numerical 


1.0 ms 


0.994 


e r f,e^ w ,A/n rf e [-.05, .05] 


Fully Numerical 


5.0 ms 


0.993 


A/U r f e [-.1, .l},e rf , e^ w = 


Fully Numerical 


5.0 ms 


0.992 



TABLE I. Summary of results. 



B. Tomographic adressing 

As discussed in the introduction, the tools of ensemble control allow us to design wave- 
forms for tomographic addressing. For example, a spatial gradient of the detuning imposed 
through an external field can be used to spatially address different regions of an ensemble. 
In our system, large gradients can be achieved by the fictitious magnetic field produced 
by the light shift associated with a circularly polarized laser beam [21]. A local unitary 
transformation can be designed as a function of the intensity of the addressing laser beam. 
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As a proof-of-principle, we consider a cold atomic gas and prepare two different "non- 
classical" states (i.e. not spin-coherent states) in two spatial regions by applying global 
control pulses in the presence of a detuning inhomogeneity. In this example we have chosen 
to synthesize = (\F = 3, nip = —3) + \F = 3, nip = 3)) / y/2 in one half of the gas and 
1^2) = \F = 3, rriF = 0) in the other. Such a distribution is highly nonclassical and in itself 
would be an interesting starting point for the study of dynamics in spinor gases [27] . In the 
ideal case, we would apply a step function to perfectly select two regions with two distinct 
detunings through, e.g., application of a light-shift phase mask. In practice, however, the 
laser field which creates the level shifts cannot be focused to a perfectly sharp edge. We 
take this into account by modeling the spatial variation in detuning by 

A(x) = A (l-e-(^ m ). (19) 

In the limit m — > 00 we recover a step function. As an example we take a = 0.5 mm and 
m — 8, which gives a transition region that is quite large compared to the wavelength of 
light and thus can be easily implemented with available optics. As depicted in Fig. Q, we 
choose Ao = 300 Hz, which is large enough that the two regions can be easily distinguished, 
but is not large compared to the rf /microwave Lamor/Rabi frequencies. 

Our goal is to synthesize in the region with detuning A = and | i n the region 
with A = 300 Hz. We have included some spread in the detuning to account for possible 
noisy background fields, and optimize on a grid around the desired detuning, A = A + 5, 
5 = {0, ±10Hz}. In this case the objective function is, 

JF(<&) = J2 (I (H U(A = 5, d>) |t) | 2 + U(A = 300Hz + 5, d>) |t) | 2 ) • (20) 

s 

In the transition region between the two targets, we do not constrain the solution. Using 
the same gradient search, we find the control waveform for the phases of the oscillating fields 
with their amplitudes fixed at the maximum values. 

The results of the optimization are shown in Fig. Q. To evaluate the performance we 
have plotted both (F z ) and (F^) as a function of position. For the two target states we 
have (Vil-FrhM = (HFzlfa) = and = 9 while (^ 2 | F 2 Z \f 2 ) = 0. In the 

region between x = and 0.3 mm, where the detuning is between A/2n = and 10 Hz, 
\(F Z )\ < 0.0125 and (F 2 2 ) > 8.95 as expected for \ipi). In the region between x = 0.6 and 
1 mm, where the detuning is between A/2ir = 290 and 300 Hz, \(F Z )\ < 1.5 x 10~ 3 and 
(F z 2 ) < 0.012, as expected for \ip 2 ). 
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(a) 300 



200 



A(Hz) 



100 




FIG. 4. Spatially selective control, (a) Plot of the spatial variation in the detuning, A, imposed by 
a light-shift gradient. The gradient separates the gas into two spatial regions in order to generate 
two distinct target states : for x < 0.5 mm, the target state is = ^(\F = 3,rnp = —3) + 
\F = 3,mp = 3)) and for x > 0.5 mm the target state is \ip2) = \F = 3, mi? = 0). Because a 
perfectly sharp transition between the two detunings is not possible in the lab, there is a transition 
region around x = 0.5 mm. The microwave and rf fields are applied simultaneously with maximum 
power and numerically optimized phases, (b) Results of the sythesized states. The left axis shows 
(F%) vs. position and the right axis shows (F z ) vs. position. If the target states are synthesized 
perfectly, then <-0i I F z \ipi) = F z \ip 2 ) = and (fa] F 2 Z \^ x ) = 9 while (ip 2 \ F% ^2) = 0. With 
the exception of the transition region at x = 0.5 mm, these goals are achieved. 

V. SUMMARY 

We have studied the control of <i-level quantum systems or qudits encoded in the hyper- 
fine magnetic sublevels of alkali-metal atoms when there are variations in the external rf and 
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microwave fields that drive the system. Such variations could be the result of uncertainties 
that result from experimental imperfection or intentionally applied in order to tomographi- 
cally address specific members of an inhomogeneous ensemble. We restricted our attention 
to qudits formed from the subspace of the lower hyperfine manifold F± = I — 1/2, plus one 
sublevel in the upper manifold Ff = 1 + 1/2; for 133 Cs this is an 8D Hilbert space. With this 
structure, we can achieve control on the system through a simple series of SU(2) rotations 
driven alternately by resonant rf and microwave fields. 

We studied the simplest control problem - state-synthesis - whereby a known fiducial 
state is transformed to an arbitrary state in the Hilbert space. Our construction is based on 
the protocol of Law an Eberly [26], extended to F > 1/2 representations of SU(2). Moreover, 
because it consists of a series of SU{2) we can draw upon the results of [HI [19], which 
proved that arbitrary unitary transformations on two-level systems could be constructed as 
a function of variations in the detuning and Rabi frequency of the driving Hamiltonian, and 
we extended them to our case of inhomogeneous control of qudits. 

We proved that the semi-analytic protocol is efficient, though not necessarily optimal in 
the duration of the pulse sequence. More importantly, this procedure provides a foundation 
for numerical searches for more efficient pulse sequences. This situation is similar to that 
found in two-level control, where a Lie algebraic approach provides a proof of principle that 
appropriate controls exist, while numerical optimization is used to find the optimum pulse 
sequence [HI [19]. The search is based on a simple gradient ascent with the objective of 
maximizing the fidelity of the target, averaged over the inhomogeneous parameters. We 
compared the performance of semi-analytical state synthesis to fully numerical state syn- 
thesis, and studied robust control over a range of experimental uncertainties. For small 
inhomogeneities, < 1%, both approaches achieve average fidelities greater than 0.99. How- 
ever, the fully numerical approach can find control sequences which require significantly 
less time. Thus, for the remainder of the paper, we studied the fully numerical approach. 
For the parameters we studied, fidelities greater than 0.99 are possible with 5% errors in 
detuning and amplitudes of the driving fields. As a testament to the efficiency of the search 
procedure, we were able to find control waveforms that are robust to 10% errors in detuning, 
well beyond experimental uncertainty. 

In addition to robust control we performed a proof-of-principle test of tomographic ad- 
dressing using designed spatial variations in the detuning over an extended ensemble. Such 
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tomography has been employed to address individual atoms in optical lattices, and con- 
sidered even under circumstances beyond the diffraction limit of an addressing laser beam 
[25]. Prior work assumed that the non-addressed members were unaffected solely because 
they were too far off resonance. Here we studied more general extensions, employing the 
tools of ensemble control. In practice, different members of the ensemble can be made 
to undergo different unitary transformations depending on the local parameters. As a 
simple example, we showed how one can synthesize, with high fidelity, the state \ipi) = 
(\F = 3, m F = — 3) + \F = 3, m F = 3))/\/2 in one half of the gas and \%jj 2 ) = \F = 3, m F = 0) 
in the other, limited by the resolution of the spatial gradient of detuning. 

In the future, we intend to extend this work in a number of directions. We have previously 
shown how one can leverage off of the efficiency of state-synthesis to design the a full unitary 
transformation of the <i-dimensional Hilbert space [Bj . For such a procedure to work, it will 
be critical take into account inhomogeneities in the driving fields. The tools of ensemble 
control developed here provide the necessary foundation. Additionally, we plan to extend 
the work of tomographic control towards addressable control of atoms in optical lattices, a 
key ingredient in many studies of quantum computing and quantum simulation. 
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